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ABSTRACT 

The magnetic field of tlie solar wind near the Sun is very difficult to measure directly. 
Measurements of Faraday rotation of linearly polarized radio sources occulted by the 
solar wind provide a unique opportunity to estimate this magnetic field, and the 
technique has been widely used in the past. However Faraday rotation is a path integral 
of the product of electron density and the projection of the magnetic field on the path. 
The electron density near the Sun can be measured by several methods, but it is quite 
variable. Here we show that it is possible to measure the path integrated electron 
density and the Faraday rotation simultaneously at 6-10 Rq using millisecond pulsars 
as the linearly polarized radio source. By analyzing the Faraday rotation measurements 
with and without the simultaneous electron density observations we show that these 
observations significantly improve the accuracy of the magnetic field estimates. 
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1 INTRODUCTION 

Tlie magnetic field in the solar corona and the inner so- 
lar wind is very important because it controls the struc- 
ture and dynamics of the coronal plasma. Unfortunately it 
is very difficult to measure the magnetic field in this re- 
gion directly. Optical observations can be made in the pho- 
tospher e, and these can be extrapolated outwards into th e 
corona (jSchatten et all Il969l : lAltschuler fc Newkirlj , 
but such extrapolations do not agree well with estimates 
made by extrapolating i nwards space-craft measurements at 
Mercury (iBurlagalbOOll '). 

It is possible to estimate the magnetic field in 
the quiet solar wind by measuring the Faraday rota- 
tion of linearly polarized radio waves propagating through 
the region of interest. Thus a good deal of work has 
been done using the Faraday rotati o n technique (e.g. 
Levy et al.l Il969i: IPatzold et al.' 'igST"; 'Bird & Edenhofei' 
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Sakurai fc Spanglcr 1994a b; Mancuso & Spanglcr 



1999 , l2000l : |jensen~et al.ll2005l : ISpangleiil2005l : llnglebv et al.i 



20071 ). The Faraday rotation of the position angle of the lin- 
ear polarization (in cgs units) is given by 



All; = y 



ne{l)B(l)-dl 

path 



(1) 



The magnetic field of the inner corona is dominated by 
closed loops over active regions and open field lines over 
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coronal holes (|Guhathakurta fc Fished 1199^ ). Closed loops 
seldom extend past 2-3 Rq and the magnetic field, as defined 
by bright striations in white light observations, appears to 
be ra dial outside of this distance (|Guhathakurta fc Fisher! 
Il995h . At larger distances a tangential component builds 
up due to rotation of the Sun (Hundhausen, 1972), but in 
the regions of interest for Faraday rotation observations the 
magnetic field is predominantly radial. This has an impor- 
tant consequence for Faraday rotation observations. If the 
corona were spherically symmetric the sign of B-dl would 
reverse at the closest point of approach and the Faraday ro- 
tation would approach zero. So one cannot use the simple 
model of a spherically symmetric solar wind - Faraday rota- 
tion is a measure of the deviation from spherical symmetry. 
A further problem with interpreting Faraday rotation obser- 
vations is that the electron density must be known. Fortu- 
nately the electron density can be measured by several meth- 
ods and its average behavior is reasonably well understood 
(e.g. lBird et al.ll 19941 : ICuhathakurta et al.lll99^ ). However it 
is time variable, and the path integration must be modeled 
carefully to obtain reliable estimates of the magnetic field. 

With a pulsar one can also measure the group delay of 
the pulse passing through the solar wind. This is given by 



A" 



2-KmeC^ 



path 



(2) 



We will show that this constraint is very valuable because 
it allows us to estimate the ratio (F) of the instantaneous 
electron density ne{t) to the average (rie). In our observa- 
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tions this factor lay in the range 0.64 < F < 1.77. Cor- 
rection of the magnetic field estimate by this factor is a 
significant improvement. However the group delay due to 
the solar wind is small and cannot be measured with suf- 
ficient accuracy for most pulsars. Fortunately the class of 
pulsars with millisecond periods (MSPs) has very high ro- 
tational stability and can be timed with a precision of the 
order of 100 ns (e.g. iManchesteil l2008l ). This is sufficien t 
to measure the solar wind contribution (|You et al.ll2007al '). 
Simultaneous measurements of Faraday rotation and group 
delay can also be performed with suitably equipped space- 
craft and such measurements have been very useful, but they 
not bee n possible since the I felios probes which worked up 
to 1985 (jHollweg et aLlligsd ). 

Coronal modeling can be done most reliably near so- 
lar minimum when the magnetic field structure is relatively 
simple and stable. At this phase of the solar activity cycle 
it can be characterized by a "magnetic pole" which may be 
displaced from the rotational pole. The magnetic field is as- 
sumed to be radial everywhere (and decreases with distance 
like to conserve flux) but is opposite in the opposite 
hemispheres. Thus the magnetic equator is a current-sheet 
where the field reverses. The current sheet appears above 
and below the rotational equator as the Sun rotates. The 
solar wind velocity is low and the electron density is high 
in a belt of about ±20° around the current sheet. However, 
the magnitude of the magnetic field is roughly the same 
in the fast and slow wind. Figure [T] shows the location of 
the current sheet and the slow wind belt during one of our 
observations. Here the current sheet, as estimated using the 
Wilcox Solar Observatory (WSO data, is shown as a heavy 
solid line and the slow wind belt is bounded by dash-dotted 
lines. Points on the line of sight from the Earth to the pul- 
sar at the time of observation are shown as a line of dots 
for which the spacing is 5° subtended at the Sun. One can 
see that the line of sight crosses the current sheet twice be- 
tween the Earth and the closest point of approach and it 
passes through two different low and high density regions. 
Thus there is a current sheet at the magnetic equator where 
the field reverses. 

However the location of the current sheet is subject 
to some error because a current-free potential-field model 
was used to extrapolate the photospheric observations of 
the WSO into the corona. Also the width of the slow wind 
belt is not well-defined and the density may vary with time. 
In this work we are concerned with the effect of these model- 
ing errors on the final magnetic field estimate, and how these 
can be reduced using the additional constraint provided by a 
si multaneous group delay measurement. This was suggested 
by lOrd et al.l (20Q7.), but they were unable to measure the 
group delay with sufficient accuracy in their observations. 
We do not have enough measurements to make a strong 
statement about the magnetic field of the solar wind, other 
than to show consistency with other measurements, but we 
have sufficient data to show the modeling errors one can 
expect and to compare them with statistical errors. In the 
following sections we will discuss: the observations and the 
primary analysis; the coronal model and the fitting process; 
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Figure 1. Tlic current siieet as determined by the WSO with 
the "classic" extrapolation to a source surface at 2.5 Rq is shown 
as a heavy solid line. The location of the slow wind belt (±20° 
about the current sheet) is indicated by dashed lines. The path to 
PSR J1022-I-1001 on 2006 Aug 25 at 00:20 UT has been extrap- 
olated back to the source surface along the streamlines assuming 
radial flow and a velocity of 400 km/s. The source-surface point 
corresponding to the point of closest approach to the Sun along 
the line of sight is indicated by a triangle and that corresponding 
to the Earth position by ®. The furthest point from the Earth 
which is shown subtends an angle of 85° at the Sun. 



the magnetic fields estimates and their sensitivity to model 
errors; and the prospects for future such observations. 



2 OBSERVATIONS AND PRIMARY ANALYSIS 

The observations were part of the Parkes Pulsar Timing Ar- 
ray (PPTA) project which commenced in 2004 and has made 
observations of 20 MSPs at intervals of 2-3 weeks since then 
(jManchesteill2008l i . Here we restrict our analysis to the pul- 
sar which provides the best observations for our purpose, 
PSR J1022±1001 which has ecliptic latitude fee = -0.06°, 
but two other MSPs regularly observed with the PPTA, 
J1730-2304 {be = 0.19°) and J1824-2452 (fee = -1.55°) 
are potentially useful. The observations analysed are those 
when the line of sight to PSR J1022±1001 happened to be 
close to the Sun - they were not deliberately scheduled to 
be near the Sun. As a result there are few observations use- 
ful for this work. All observations us e the center beam of 
the P arkes 20-cm Multibeam receiver (|Stavelev-Smith et al.l 
Il996t ) at frequencies close to 1400 MHz. All of the data 
were recorded in full-polarization mode. The back-end sig- 
nal processing systems are Parkes digital filter bank sys- 
tems (PDFBl to PDFB4). Each observation has a total 
bandwidth of 256 MHz with 1024 channels. The duration of 
each observation is 64 minutes and the data were averaged 
with 1-min sub-integrations. A 2-min pulsed calibration sig- 
nal was recorded before each pulsar observation. The flux 
density scale was set using observations of Hydra A. The 
cross-coupling between feed probes was measured using ob- 
servations of PSR J0437— 4715 covering a wide range of hour 
angles. 

The primary analysis was done using P SRC HIVE soft- 
ware (H otan et al. 2004) . Initially, 5% of the band edges and 
radio frequency interference (RFI) were automatically re- 
moved. Then the data were flux- and polarization-calibrated 
and the cross-coupling between the feed polarizations was 
corrected using the pac/pcm algorithm in PSRCHIVE. In 
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pulsar analysis the group delay is quantified by the "dis- 
persion measure" DM which is the path integrated electron 
density in pc cm~^. The Faraday rotation is quantified by 
the "rotation measure" RM, which is the term in square 
brackets in equation 1 in units of rad m~^. We will use these 
measures in our further discussion. 

The method of obtai ning RM is similar to that dis- 
cussed in lYan et all l|20lj ). We approached the final result 
with a series of three successive approximations. We first ad- 
justed the RM to optimize the linearly polarized intensity 
integrated over the bandwidth. Then we separated each ob- 
servation into upper and lower halves of the total band. An 
improved value of RM was determined using the weighted 
mean of position angle difference between the two halves of 
the band for each pulse phase bin. The RM also includes a 
time variable ionospheric component. We corrected it using 
an ionospheric model which is encoded in 2007 International 
Reference lonophere (IRI) modeQ. It was relatively small, 
ranging from —0.45 to — 2.25radm~^, but not negligible. 

To obtain the final approximation to RM we used the 
corrected p osition ang l e vs p ulse phase measured far from 
the Sun by lYan et aU (|201ll ') as a reference template. We 
then found a weighted average of the difTerences of each 
observation near the Sun to this reference template. The 
normalized chi-squared of these weighted averages was typ- 
ically about 1.3, indicating that the error estimates used in 
the weighted fit were reasonably accurate and the position- 
angle template was a good model when the pulsar was near 
the Sun. 

To estimate the DM contribution from the solar wind, 
we used a simpler procedure. The DM variations due to the 
interstellar medium for PSR J1022+100 1 are smaller than 
the measurement error (IYou et al.ll2007bl ). So we simply es- 
timated pulse arrival times in the normal way, by fitting 
a detailed timing model to the entire set of observations, 
excluding those near the Sun. We then obtained DM© by 
comparing the observed timing residual near the Sun with 
this distant average. There is a contribution from the iono- 
sphere but, unlike the RM contribution, the ionospheric DM 
contribution is negligible compared with the measurement 
error. 

The final RM and DM measurements used are summa- 
rized in Tableland plotted in Figure[2]vs distance from the 
Sun. 



3 SOLAR WIND MODEL FITTING 

We assume that the electron density is bimodal, as 
noted earlier. The average density in the two modes 
has been estimated ([G uhathakurta fc F isher! 1 19951 . 1 19981 : 
iMuhleman fc Anderson.. 1981 : Allen 1947i) ; we us e the same 
approximations as were used bv lYou et al.l l|2007al 'l . The den- 
sity rie in the fast wind is 

and Tie in the slow wind is 

UeiR) = {4.1R-^ + 23.53i?-^'^)10" + 

(1.5/?-^ + 2.99J?-'^)10'^ m-^ (4) 

^ See |http: / /lrl.gsfc.nasa.gov/| 
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Figure 2. Observations of PSR J1022+1001 near the Sun. The 
abscissa is the distance from the Sun. The E limb data is plotted 
with a negative distance, and the W limb data with a positive 
distance. The observations for different years are marked by dif- 
ferent symbols. The error bars on the RM are smaller than the 
symbol size. 

Here R is the distance from the center of the Sun in _R0 . To 
match the DM observations we scale the electron density by 
a factor F, and to minimize the number of free parameters 
we use the same factor in both the fast and slow wind. 

We assume that the slow wind belt is centered on the 
current sheet as determined by the W SO pho tosphcric obser- 
vations extrapolated into the corona (|McComas et al. 2000|). 
However the width of the slow wind belt is variable and 
not precisely determined in any case (|Mancuso fc Spangled 
^2000). Accordingly we have tested models with widths of 
±10°, ±20° and ±30°. 

The extrapolation is done with three similar, but subtly 
different, techniques referred to as: "classic"; "radial 250"; 
and "radial 325" in the WSO data bascEI. The three models 
find a potential solution which is current free and radial at 
some source surface. The radial models both constrain the 
photospheric field to be radial and differ only in the source 
surface which is at 2.50 or 3.25 Rq. The classic model does 
not constrain the photospheric field to be radial but requires 
a somewhat ad hoc polar field correction. It assumes a source 
surface at 2.5 Rq. We assume that the magnetic field is the 
same in the fast and slow wind, that it is radial, and that it 
varies quadratically with solar distance. We take the polarity 
of the magnetic field from the corresponding WSO model. 

To provide an idea of the range of the extrapolation 
models and a comparison with other solar wind observations, 
we plot the current sheets for the three models over the solar 
wind velocity measurements for 2006 from the STELAlfl 
in Figure [S] Here one can see that the three extrapolation 
models vary significantly in the maximum latitude extent of 
the current sheet, and that the width of the slow velocity 
belt is not clearly defined. 

It is helpful to examine the integrand of the path in- 
tegrals to see where the maximum contributions arise. For 
numerical purposes we transform the path integral into an 

^ see |http: / / wso.stanford.edu/synsourcel. html| 

bee |http:// www. stelab.nagoya-u. ac.jp/omosaic/crle. htmll 
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Table 1. Observations of the PSR J1022+1001 near the Sun. The echptic coordinates are be = 



-0.06° and le = 153.9°. 



Date 


UT 


Ro 


/ 


L/I 


RMq 


DMq 






(R©) 


(mjy) 


(%) 


(rad m~-^) 


(10~^cm~'^pc) 


2005/08/29 


02:15 


7.5 


3.0 ±0.3 


62 


-11.6 ±0.15 


3.6 ± 0.8 


2006/08/24 


03:40 


-11.3 


8.0 ±0.4 


52 


-0.72 ±0.06 


1.2 ±0.3 


2006/08/25 


00:20 


-8.2 


4.6 ±1.0 


48 


-1.49 ±0.24 


0.3 ±2.0 


2006/08/28 


23:49 


6.2 


13.3 ±1.3 


58 


-23.3 ±0.09 


8.2 ±0.6 


2006/08/30 


00:02 


9.9 


12.4 ±0.7 


52 


-3.9 ±0.05 


3.6 ±0.3 


2006/08/31 


01:19 


13.7 


2.3 ±0.4 


54 


0.34 ±0.16 


2.0 ± 1.3 


2008/08/31 


23:41 


19.0 


1.1 ±0.2 


48 


-0.18 ±0.17 


5.6 ±0.9 


2009/08/30 


02:01 


11.2 


3.2 ±0.2 


55 


0.17 ±0.08 


7.7 ±0.4 
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Figure 3. Solar wind velocity from STELAB, Nagoya U. This 
figure has been rotated in longitude to match Figure [T] The cur- 
rent sheet locations from the three WSO models have been over- 
laid. The "classic" model is a solid line, the "radial 250" model 
is dashed, and the "radial 325" model is dot-dashed. 
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Figure 4. Integrands for group delay (top panel) and Faraday 
rotation (bottom panel) for the path shovirn in Figure [T] The 
integrands for fast and slow wind alone are shown dashed in the 
top panel. The integrand for Faraday rotation if the phase reversal 
at the current sheet were not present is shown dashed in the lower 
panel. The X marks the point near Carrington longitude 310° 
where the line of sight crosses the current sheet (cf. Figure [TJ. 



integral over the angle subtended at the Sun. This makes 
the range finite and compresses the integrand in a useful 
way. The integrands for group delay and Faraday rotation 
for the path shown in Figure [1] are shown in Figure ID One 
can see that the contributions from the slow dense wind are 
dominant in both integrands and the location of the current 
sheet is very important in obtaining the Faraday rotation. 



The fitting procedure consists of three steps. First, we 
examine the LASCO C3 movie^on the three days surround- 
ing the observations and confirm that there was no obvious 
transient during the observing period. This was true of all 
the observations in Table [T] Second, we calculate the DM 
from one of the 9 models (3 extrapolation methods and 3 
widths of the slow wind belt) and find the factor F necessary 
to match the observed DM in each case. Third, we use that 
factor F and an arbitrary \B\ to calculate the RM for that 
model, and scale \B\ until the model matches the observed 
RM. The DM and RM steps are repeated for all 9 models. 

There were three observations that we could not model 
completely in this way. In one case, 2006 Aug 25, the DM 
measurement is not adequate to determine the factor F, so 
we set F=l. In the second case, 2009 Aug 30, all three models 
showed that the path was entirely in the slow wind. The 
electron density was 20% above its average (F = 1-2), but 
the RM was very small - within 2 ct of zero. This would be 
expected if the path did not cross the current sheet. Two 
of the three extrapolation models indicated that the path 
would cross the current sheet and give a significant RM, but 
the third, the "radial 250" model did not indicate a current 
sheet crossing so the predicted RM = 0. In this case we have 
a consistent model but no estimate of the magnetic field. In 
the third case, 2008 Aug 31, all three models show that 
the path is almost entirely in the slow wind. The electron 
density was high (F = 3.4 ±0.4) and the RM was very small. 
The "radial 325" model predicted the wrong sign of RM 
indicating a significant error in the location of the current 
sheet. The "classic" model predicted \B\ — 0.2 ±0.2 mG and 
the "radial 250" model |_B| = 59 ± 56 mG. We doubt that 
\B\ < 0.4 mG at 11.2 Rq and suggest that the "radial 250" 
model is more consistent with the observation. 

The results for the nine different models are shown in 
Table [2] One can see that the factor F varies from 0.5 to 
3.0, which confirms that having simultaneous DM observa- 
tions provides a strong constraint on the solar wind model. 
We have taken the mean of the 9 estimates of for each 
observation as the best estimate of for a given observa- 
tion. The rms of these 9 estimates is a measure of the model 
sensitivity of that particular estimate. We also have an esti- 
mate of the statistical error derived from propagation of the 
errors on the primary measurements of DM and RM. These 
statistical errors are dominated by the error in measuring 
DM. Finally we find the \B\ estimate that would have been 

^ Seejhttp://soho.esac.esa.int/data/archive/index_ssa.html| 
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Table 2. Results of model sensitivity analysis using 3 WSO mod- 
els and 3 widths for slow wind. Here \B\ in mG is estimated at 
the point where the line of sight is closest to the Sun. 



Observation 


mod 


±10° 


±20° 


±30° 






F 


\B\ 


F 




F 


\B\ 


2005/08/29 


clas 


0.9 


48.9 


0.6 


33.1 


0.5 


26.6 


02:16 UT 


250 


0.9 


29.6 


0.6 


23.2 


0.5 


20.0 


7.5 Rq 


325 


0.8 


23.7 


0.5 


19.3 


0.5 


19.0 


2006/08/28 


clas 


3.0 


24.2 


2.1 


21.9 


1.3 


19.8 


23:49 UT 


250 


2.7 


21.8 


1.4 


17.2 


0.9 


14.8 


6.2 Rq 


325 


2.4 


20.8 


1.3 


16.4 


0.8 


15.0 


2006/08/30 


clas 


2.3 


12.0 


1.4 


12.7 


1.0 


10.4 


00:02 UT 


250 


1.7 


10.2 


1.0 


8.2 


0.7 


6.5 


9.9 Rq 


325 


1.8 


8.4 


0.9 


7.4 


0.6 


6.7 


2006/08/25 


clas 


1.0 


14.1 


1.0 


3.5 


1.0 


2.7 


20:56 UT 


250 


1.0 


28.5 


1.0 


3.4 


1.0 


3.4 


-8.2i?o 


325 


1.0 


10.6 


1.0 


2.3 


1.0 


1.2 



made in the absence of DM observations simply by setting 
F = 1 for each model and recalculating \B\. This summary 
is shown in Table |3] Here one can see that the model errors 
(merr) exceed the statistical (sterr) errors in every case, and 
the average effect of correcting for the observed DM sub- 
stantially exceeds the model errors. 



4 DISCUSSION OF RESULTS 

We have identified four primary contributions to the error 
on magnetic field estimates made using the Faraday rota- 
tion technique. In order of importance these are: failure to 
correct for DM variations (~50%); inaccuracy of estimation 
of the current sheet location («25%); error in the DM mea- 
surement («15%); and error in the Faraday rotation esti- 
mate Thus use of simultaneous DM observations can 
roughly double the precision of the magnetic field estimates. 
The value of the PPTA for such measurements is clear. It 
provides a good timing model for the DM far from the Sun 
and a good template for the position-angle variations far 
from the Sun. 

We have plotted our estimates of along with those of 
oth er authors in FigureEl T he estimates of iOrd et al. ( 2007 ) 
andU mirnova et al.1 12009) are not plotted because in both 
cases the authors only provide a measurement of the average 
line-of-sight component of the magnetic field and do not at- 
tempt to estimate \B\. Their estimates of the average line-of- 
sight components are consistent with other observations but 
do not put a useful constraint on \B\. Our estimates of |_B 
are somewhat lower than the others, but are consistent with 
them, within the error bars. All the observations are broadly 
consistent with the Helios space craft measurements and an 
extrapolation like R~^. Thus our measurements, although 
few in number, confirm the very important result that the 
potential field models of the corona seriously underestimate 
the magnitude of the magnetic field for R > 5Rq. 

We did not observe any transients near the Sun and 
would not have been able to analyze them with the tech- 
nique we used. However it might be possible to do so if 
sufficient additional information, such as might be provided 
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Figure 5. Magnetic field observations near the Sun. The mea- 
surements in Table \3\ are plotted as * symbols with er ror bars. 
The extrapolated Helio s measurements are a solid line l|BurIagal 
2001), an average of the lPatzold et al.l lll987l ) observations {\B\ = 
1000(6iJ-3 + I.ISR-^) mG) are a dot-das hed line and the ob- 
servations of iGopalswamv fc Yashirol ll201ll'l are + symbol s con- 
nected by a dashed line. The obs ervation of | Sp anp|ler |20oJ is an 



open box and the range given bv llnglebv et al. is a pair of 

connected x symbols. 

by STEREO were available to constrain the geometry of 
the magnetic field. We note that Faraday rotation is capable 
of 1% precision by itself, and one might be able to measure 
rapid variations of magnetic field during transients with this 
accuracy. 

We believe that more extensive observations of PSRs 
J1022+1001, J1730-2304, and J1824-2452 near the Sun, 
perhaps involving more telescopes to provide continuous 
coverage during the closest approach, would be valuable and 
should be undertaken at the next opportunity. There are 
five additional MSPs with ecliptic latitudes < 2° that may 
be suitable for precision timing: J0030+0451; J1614-2230; 
J1721-2457; J1802-2124; and J1811-24. Preliminary ob- 
servations would be required to confirm that these are "good 
timers" before undertaking coronal measurements. The clus- 
ter of pulsars with a right ascension close to 18:00 would fa- 
cilitate an intensive observing program in the last two weeks 
of December each year, but there are insufficient MSPs to 
monitor the corona on a regular basis. 
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Table 3. Final estimates of |B| at the point where the line of sight is closest to the Sun. |S| is given both corrected by measured DM 
and uncorrected. The measured DM exceeds that expected from the electron density model by the factor F. 



Observation 


Ro 




\B\ (m 


,G) 






F 






(Re) 


DMcor 


uncor 


merr 


sterr 


mean 


merr 


sterr 


2005/08/29 


7.5 


27.0 


42.3 


±9.5 


±6.7 


0.64 


±0.17 


±0.14 


2006/08/28 


6.2 


19.1 


12.3 


±3.4 


±1.4 


1.77 


±0.80 


±0.13 


2006/08/30 


9.9 


9.2 


8.0 


±2.3 


±1.1 


1.27 


±0.57 


±0.11 


2006/08/25 


-8.2 




7.7 


±8.9 


±1.5 


1.0 
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